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Computed  properties  are  presented  that  characterize  electron  emission  with  a termi- 
nated  image  potential  for  fields  from  10°  to  10°  volts  per  centimeter,  temperatures  from 
0*^  to  3000*^  K,  work  functions  from  1 to  8 volts,  and  Fermi  levels  from  1 to  1 5 electron 
volts.  Theoretic  electron- number  and  current  densities,  energy  distributions  and  fluxes, 
and  temperatures  describe  the  electrons  that  escape  over  and  through  confining  potential 
walls.  Because  the  terminated  image  potential  produces  higher  and  wider  emission  bar- 
riers, its  current  densities  fall  below  those  for  the  ordinary  (nonterminated)  image  po- 
tential at  identical  conditions.  Where  they  apply,  these  theories  probably  bracket  real 
thermal,  field  emission. 


INTRODUCTION 

The  calculations  of  this  paper  predict  properties  of  thermal,  field  emission  impeded 
by  a barrier  with  an  image  potential  that  terminates  at  the  Fermi  level  on  the  surface  of 
the  emitter  (TIP)  (ref.  1).  These  theoretic  electron- number  and  current  densities, 
energy  distributions  and  fluxes,  temperatures,  and  barrier  dimensions  reveal  supra- 
barrier  and  intrabarrier  emission  effects.  To  reach  regions  of  thermal,  field  emission 
where  simple  theories  probably  apply  as  well  as  where  they  are  generally  used,  the  re- 
sults spread  over  fields  from  10®  to  10®  volts  per  centimeter,  temperatures  from  0°  to 
3000®  K,  work  functions  from  1 to  8 volts,  and  Fermi  levels  from  1 to  15  electron  volts. 

Dyke  and  Dolan  (ref.  2)  published  some  similar  numbers  for  the  usual  nonterminated 
image  potential  (NIP).  They  sketched  energy  distributions  for  a 4.  5-volt  work  function 
and  tabulated  current  densities  for  work  functions  of  4,  4.  5,  and  5 volts  with  fields  from 
lo"^  to  10®  volts  per  centimeter  and  temperatures  from  0®  to  3500®  K. 


The  NIP  and  TIP  begin  and  end  a series  of  simple  models  for  emission  barriers  that 
should  include  an  approximation  of  reality  (ref.  1).  Therefore,  actual  thermal,  field 
electronics  probably  interpolates  the  NIP  and  TIP  theories  where  they  hold. 

The  TIP  properties  presented  here  came  from  mathematic  processes  that  were 
transcribed  in  FORTRAN  by  Susan  Button  and  Annie  Easley  and  are  recorded  in  appen- 
dix C. 


THEORY 


The  TIP  theory  for  thermal,  field  emission  describes  the  escape  of  electrons  over 
and  through  the  barrier  profiled  in  figure  1.  Its  rationale  resides  in  reference  1,  which 
details  the  TIP  derivation.  This  section  merely  summarizes  that  development;  appen- 
dix A defines  the  symbols. 

Electrons  top  the  TIP  barrier  in  the  following  current  densities: 
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Here  the  zero-order  term  approximates  the  entire  expression  adequately  for  usual  work 
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functions  to  fields  greater  than  10  volts  per  centimeter,  where  tunneling  prevails  and 
the  theory  fails  (ref.  1) . 

Penetration  probabilities  for  the  TIP  barrier  start  with  the  WKB  expression  and  re- 
strictions (ref.  3), 


(2) 
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where  f(V,e)  approximates  unity  and  varies  slowly  (ref.  2).  Reference  1 gives  two  re- 
sults for  equation  (2),  one  each  for  the  ranges  of  barrier  potentials  above  and  below  the 
Fermi  level.  In  these  solutions  the  exponents  involve  elliptic  integrals  of  the  first  and 
second  kinds  and  yield  values  like  those  tabulated  in  reference  1. 

When  TIP  penetration  probabilities  (eq.  (2))  multiply  the  Sommerfield-Bethe  electron 
supply  function  (ref.  4), 


n(v^)dv^  = 

i 

the  product  predicts  tunneling  distributions.  This  intrabarrier  effect  and  the  suprabarrier 
equation  (1)  compose  the  total  TIP  thermal,  field  emission. 

Properties  of  the  freed  electrons  depend  on  emitter  characteristics  {ji,  <p,  and  T) 
and  the  field  (E).  Average  conditions  for  emitted  electrons  ((  u))  derive  from  their 
velocity  distribution  function,  as  the  following  expression  reveals: 


(U)  = 


{ P(Vx)n(Vx)u(Vjj)dVj 


(4) 


This  equation  integrates  both  intrabarrier  and  suprabarrier  effects  with  a penetration 
probability  [P(v^)J  dictated  by  equation  (2)  for  potentials  up  to  the  summit  of  the  barrier 
and  fixed  at  unity  for  energies  above  that  point. 


DEFINITIONS  FOR  COMPUTED  RESULTS 

When  the  information  of  the  previous  section  was  detailed,  encoded,  and  relayed  to 
the  IBM  7094,  it  poured  forth  results  from  equations  (2)  to  (4).  In  its  present  state  of 
learning,  the  machine  integrates  numerically  the  basic  penetration  expression  (eq.  (2)) 
more  readily  than  it  stores,  interpolates,  and  dispenses  data  from  tabulations  of  values 
of  a function  of  several  variables.  Thus,  the  final  solutions  for  equation  (2)  and  the  nu- 
meric results  from  reference  1 were  not  used  in  these  calculations. 

With  the  FORTRAN  statement  for  TIP  calculations  (appendix  C),  the  IBM  7094  pro- 
duced output  sheets  typified  by  figure  2 and  described  in  appendix  B.  These  results 
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comprise  emitter  conditions  (first  row);  discrete,  fractional,  and  differential  properties 
for  emission  (columns);  average  characteristics  (last  two  rows);  and,  finally,  four  plots 
of  energy  distributions  of  number  and  current  densities  (the  seventh  and  eighth  columns 
against  mean  values  of  the  first  column,  and  the  ninth  and  tenth  columns  against  the  sec- 
ond column).  In  the  general  presentation  of  results  (fig.  3)  only  the  three  rows  and  four 
plots  appear,  but  the  graphs  also  display  tabulations  of  values  for  abscissae  and  ordi- 
nates, which  for  all  four  plots  include  the  first  two  and  last  four  columns. 

For  these  calculations,  escaping  electrons  begin  in  the  emitter  and  end  in  free  space 
at  the  same  potential  (fig.  1).  Therefore,  the  tabulated  electron  properties  refer  to  the 
bottom  of  the  conduction  band;  this  potential  prohibits  no  emission  and  exhibits  no  over- 
all acceleration.  Furthermore,  the  results  describe  either  those  internal  electrons 
about  to  leave  or  the  emitted  ones. 

DISCUSSION  OF  RESULTS 

What  conditions  best  test  the  NIP  and  TIP  models?  And  where  should  these  simple 
theories  be  used? 

A high  electric  field  separates  and  emphasizes  the  NIP  and  TIP  effects;  they  diverge 
as  tunneling  increases,  while  real  thermal,  field  emission  probably  falls  between  them 
in  the  region  where  image  potentials  apply.  Apparently  then,  high  fields  promise  a good 
evaluation  of  the  NIP  and  TIP  approximations.  In  fact,  most  studies  compare  experi- 
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ments  with  theory  above  10  volts  per  centimeter. 
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At  10  volts  per  centimeter,  however,  NIP  and  TIP  barriers  maximize  about  5 to 

O 

6 angstroms  from  the  emitter  face,  and  at  10  volts  per  centimeter  these  potentials  peak 
within  2 angstroms  from  the  surface.  In  the  light  of  these  dimensions  the  assumption  of 
the  smooth,  planar  emitter  surface  stands  stark.  Where  the  thickness  of  the  potential 
wall  nears  the  size  of  superficial  imperfections,  interfacial  fields,  which  definitely  differ 
from  those  for  the  assumed  flat  face,  strongly  influence  the  overall  shape  of  the  barrier. 

Furthermore,  high  fields  probably  penetrate  the  emitter  and  pull  down  the  bottom  of 
the  conduction  band  near  the  boundary  (ref.  5).  This  potential  incline  draws  more  elec- 
trons toward  the  surface,  where  they  escape  over  and  through  the  barrier  in  greater  num- 
bers. 

In  addition,  the  large  currents  and  steep  gradients  of  conditions  near  the  interface 
destroy  the  approximate  equilibrium.  Then  the  Fermi  distribution  flounders  in  its  repre- 
sentation of  emitter  electrons. 

When  a high  field  mobilizes  these  phenomena,  they  militate  against  the  reign  of  the 
simple  emission  models.  But  when  fields  reduce  to  10  volts  per  centimeter,  the  NIP  and 
TIP  maximums  move  18  to  19  angstroms  from  the  emitter;  this  distance  is  several  times 
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TABLE  I.  - THEORETIC  CURRENT  DENSITIES  FOR  THERMAL  FIELD 


EMISSION  WITH  4-VOLT  WORK  FUNCTION 


Source 

Reference  2 

Present  paper 

Theory 

Fowler- 

Nordheim 

solution 

(FN) 

Ordinary  image 
potential 
(NIP) 

Terminated  image 
potential 
(TIP) 

Temperature, 
T,  °K 

0 

1000 

3000 

. . j 

0 

1000 

3000 

Field,  V/cm 

2 

Current  density,  A/cm 

10^-5 

4.7X10^ 

6.  5X10® 

2.6x10® 

5.  8x10^ 

l.lxlO® 

8.2X10® 

00 

O 

2.  5x10® 

1.9X10® 

2.  9X10® 

2.1x10® 

2.  3x10® 

4.2x10® 

the  atomic  radius,  which  means  that  high-field  distortions  practically  disappear.  Ob- 
viously then,  lower  fields  favor  the  NIP  and  TIP  theories,  but  in  this  region  these  simple 
models  coalesce. 

This  conflict  rationalizes  the  range  of  results  (fig.  3),  which  span  from  10®  volts  per 
centimeter,  where  field  emission  hardly  appears,  to  10®  volts  per  centimeter,  where  the 
TIP  model  hardly  applies.  Fermi  levels  from  1 to  1 5 electron  volts  and  work  functions 
from  1 to  8 volts  cover  clean  and  coated  metal  emitters,  while  temperatures  from  0°  to 

3000°  K run  the  gamut  of  practical  emission  from  solids. 

7 8 

Between  10  and  10  volts  per  centimeter,  the  TIP  results  correspond  to  those  for 
the  NIP  and  Fowler-Nordheim  (FN)  solutions  tabulated  in  reference  2.  Table  I compares 
current  densities  obtained  with  FN  (0°  K),  NIP  (1000°  and  3000°  K),  and  TIP  (0°,  1000°, 
and  3000°  K)  for  a 4- volt  work  function  and  fields  of  10^'  ® and  10®  volts  per  centimeter. 
The  TIP  values  fall  considerably  below  those  for  FN  and  NIP.  Unfortunately  in  refer- 
ence 2 at  10  volts  per  centimeter,  the  FN  value  exceeds  the  NIP  number  for  1000°  K at 
10  volts  per  centimeter,  which  precludes  a valid  comparison  for  0°  K. 

Table  I specifies  no  Fermi  level  because  the  emitter  potential  appears  not  to  affect 
current  density.  Technically  this  lacks  some  part  of  the  truth,  as  equation  (5)  indicates: 
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Except  for  high  fields  and  low  work  functions,  however,  plunging  penetration  probabilities 
choke  off  emission  at  the  low-energy  end  of  the  distribution  before  the  outward  directed 
velocities  reach  zero  (e^-  fji  = -ji)  and  cut  off  the  electron  supply.  Therefore,  equa- 
tion (3)  of  reference  2,  where  -°°  rather  than  -[i  limits  the  integration,  generally 
approximates  equation  (5)  well. 


r°° 

j«/  V{cp,E,e  -ii)de^  = K(P,T,E)  (6) 

OO 

But  the  Fermi  level  influences  other  properties  of  thermal,  field  emission  presented  in 
figures  2 and  3. 

In  reference  2 another  high- field  complication  arises  when  NIP  potential  maximums 

O 

sink  below  the  Fermi  level  near  10  volts  per  centimeter.  Of  course,  the  top  of  the  TIP 
barrier  cannot  slip  below  the  Fermi  level  because  the  image  and  emitter  potentials  join 
at  the  interface.  There  is  little  need  to  belabor  this  point  further,  however,  because 
fields  near  10^  volts  per  centimeter  overwhelm  the  NIP  and  TIP  models  at  every  turn. 

There  remains  one  apparent  anomaly  to  explain.  For  the  combination  of  low  Fermi 
levels  and  high  fields,  distribution  functions  for  emitted  electrons  display  a definite  dis- 
continuity at  e = 0.  Although  this  cutoff  always  occurs,  as  equation  (5)  reveals,  emis- 

A 

sion  usually  fades  to  negligibility  well  above  the  barrier  bottom.  In  previous  papers 
(i.  e. , ref.  2)  no  such  breaks  in  emission  occur.  Those  studies  refer  electron  energies 
within  the  emitter  to  the  Fermi  level  and  extend  the  supply  function  to  plus  and  minus 
infinity  in  an  approximation  to  ease  integration  (eq.  (6)).  That  approach  requires  no 
Fermi  level,  gives  no  energy  flux,  and  exhibits  no  discontinuity  in  emission.  In  cont- 
trast,  the  present  paper  specifies  the  Fermi  level,  yields  complete  information  on  ener- 
gies of  emitted  electrons,  and  cuts  off  emission  at  = 0.  Appreciable  discontinuities 
in  TIP  emission  appear,  however,  only  at  fields  that  badly  distort  the  model.  No  demon- 
strable theoretic  difficulties  appear  at  conditions  suitable  for  applications  of  NIP’s  and 
TIP’s. 

But  some  problems  arise  because  the  calculations  and  plots  are  incremental  rather 
than  continuous.  First,  where  significant  results  occur  in  regions  having  Av^  similar 
in  size  to  v^^  (at  E = 10®  V/cm  and  m = 1 eV),  the  fidelity  of  the  computing  procedure 
fails.  Again  though,  this  distortion  develops  only  at  conditions  that  preclude  the  use  of  TIP 
and  NIP  models.  Second,  although  the  distribution  functions  should  yield  smooth  curves, 
some  of  the  graphs  look  lumpy.  This  results  from  inadvertent  correlations  of  computing 
and  graphing  increments  that  prevent  points  from  scattering  randomly  about  the  true  mean 
tendencies  of  the  plotted  functions.  Thus,  a few  curves  shift  slightly  or  bend  abruptly 
where  the  distributions  actually  continue  smoothly.  In  any  event,  tabulations  on  each  of 
figures  3 indicate  this  graphic  difficulty  and  supply  proper  values  of  the  variables. 
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After  these  apologies,  the  results  confer  some  conveniences.  First,  the  NIP  and  TIP 
theories  stand  as  about  the  only  two  choices  to  predict  thermal,  field  emission;  the  value 
of  this  observation  depends  on  the  viewpoint.  Second,  realistic  situations  occur  where 
neither  penetration  currents  nor  the  models  fail;  here  the  theories  submit  tractably  to  ex- 
perimental tests.  Third,  the  computed  distributions  and  average  results  provide  bases 
for  fractional;  differential;  or  gross  evaluations  of  thermal,  field  emission,  using  re- 
tardation and  calorimetry.  Fourth,  the  spread  of  variables  allows  interpolation,  and  the 
FORTRAN  statement  (appendix  C)  facilitates  calculation  of  results  for  specific  conditions. 


CONCLUDING  REMARKS 


Although  this  report  records  results  for  fields  from  10®  to  10®  volts  per  centimeter, 
it  cautions  that  models  for  thermal,  field  emission  opposed  by  barriers  based  solely  on 
image  potentials  hold  little  hope  for  applicability  near  the  top  of  this  range.  While  the 
NIP  and  TIP  theories  may  perform  poorly  between  lo"^  and  10®  volts  per  centimeter,  they 
should  prove  adequate  up  to  10®‘  ® volts  per  centimeter.  This  assumption  gains  stature 
from  the  good  performance  of  Schottky’s  emission  equation  at  fields  near  10®  volts  per 
centimeter. 

Therefore,  experiments  to  test  and  attempts  to  apply  these  simple  models  of  emitters 
should  be  confined  to  the  low- field  region  of  thermal,  field  emission.  Refined  theories 
recognizing  the  surface  problems  that  dominate  the  potential  barrier  at  high  fields  must 
be  developed  to  predict  thermal  field  emission  above  lo"^  volts  per  centimeter. 
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APPENDIX  A 


SYMBOLS 


Equations  have  cgs  units;  results  are  given  in  units  likely  to  be  used  in  experiments. 


E 

e 

f(  ) 
h,K 

j 

m 

NIP 


n(€^),  n(€ 
n(v^) 


To 

TIP 

{^) 

u(v^) 

V 


X 


electrostatic  field 
electron  charge 

function  of  parenthasized  variables 
Planck’s  constant  and  h divided  by  2n 
current  density 
mass  of  electron 

nonterminated  (ordinary)  image  potential,  also  emission  model,  theory, 
and  other  items  related  to  or  based  on  NIP 

number  densities  of  emitted  and  emitter  electrons 

emitted-electron  number  density  distributions  as  functions  of  total 
kinetic  energy,  kinetic  energy  based  on  outward  directed  velocity 
components  only,  outward  directed  velocity  component 

penetration  probability 

emitter  temperature 

absolute  static  temperature  of  an  emitted  electron  beam  that  maintains 
its  average  velocity  and  randomizes  the  other  velocity  components 
around  its  drift  velocity 

absolute  total  temperature  of  emitted  electron  beam 

terminated  image  potential,  also  emission  model,  theory,  and  other 
items  related  to  or  based  on  TIP 

average  value  of  u 

a function  of  v^^ 

potential  relative  to  bottom  of  conduction  band 
positive  x-directed  velocity  component 
dummy  variable  for  summation  over  v^^ 

direction  normal  to  emitter  face,  positive  from  surface  into  free  space 
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inner  and  outer  turniii^  points  on  emitter  potential  barrier  (where  eV  - €^^=0) 
M + (e<p  - €^) 

total  kinetic  energy  of  an  electron 

electron  kinetic  energy  based  only  on  positive  x-directed  velocity  component 

electron  kinetic  energy  based  on  y and  z components  (parallel  with  emitter 
face)  of  velocity 

Boltzmann  constant 

emitter  Fermi  level 

emitter  work  function 


APPENDIX  B 


First  Row 

T 

E 

PHI 

AMU 

EVMAX 

Columns 

EPSX 

EPST(VXM) 


SYMBOLS  FOR  IBM  OUTPUT  SHEETS 

emitter  temperature,  °K 
electric  field,  V/ cm 
work  function,  <p,  V 
Fermi  level,  p , eV 

electron- energy  equivalent  of  top  of  TIP  barrier,  = ix  + ecp 

- (e^E)^/^  + e^E/(4<p),  eV 


highest  value  in  each  increment  of  electron  kinetic  energy  based  only  on 
outward  (positive-x)  directed  velocity  component,  = niv^2,  eV 
The  remaining  columnar  entries  pertain  to  the  mean  value  for 
each  increment.  The  form  F(VXM)  denotes  a function  of  v^  computed 
for  this  incremental  mean  € . 

average  total  kinetic  energy  of  an  internal  electron  with 


£ 


X’ 


eV: 


£t 


= €, 


+ (eyz> 


X 


le  =0 


Cy^Ml  + e 


kT 


/ 


+ I Inll  + e 


kT 


d€„ 


yz  I In  l + e 


kT 


yz 


In  ll  + e 


kT 


In  1 + e 


kT  ^ 
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when  (e^  - /i)/kT  grows  large, 


Xi,X2 


P(VXM) 

PNDV 


e 


yz 


(^)  - 

X 


= kT 


kT 


and  when  T approached  zero, 


<^y,> 


/ 


^yz  ^x  ■ ^ 


kt 


M - e. 


^ - s 

xT 


inner  and  outer  turning  points  (at  eV  - = 0)  of  TIP  barrier,  Xj 

and  X2,  cm 


2eE  \ 


TIP  penetration  probability,  P(v^),  eq.  (2) 

fraction  of  emitted  (or  about  to  be  emitted)  electrons  that  have  positive 
x-directed  velocity  components  up  to  v 


£ P(v’ )n(v;^)Av’ 
v'  =0  ^ XX 

X 


£ P(v’ )n(v*  )Av’ 

yt  — Q A A X 


11 


PN(EX) 


J(EX)/DEX 


PN(ET) 


where  eq.  (3)  defines  n(Vj^).  The  summations  proceed  through  pre- 
scribed increments. 

distribution  of  emitted  electrons  as  function  of  kinetic  energy  based  on 
positive  x-directed  velocity  component  only,  electrons/ cm  / eV 

Pn(€x)  = [P  of  eq.  (2)][n(€j^)  of  eq.  (3)] 


rate  of  change  of  emitted  current  density  with  kinetic  energy  based  on 
positive  x-directed  velocity  component,  A/cm  / eV 


^ = ^n{e^)v^e  = [pN(EX)](VXM)e 


distribution  of  emitted  electrons  as  ftmction  of  total  kinetic  energy, 

O 

electrons/ cm  /eV 


Pn(€iji)  = Pn(Cjj) 


where 
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J(ET)/DET 


As  (m  - €^)/kT  grows  large,  (€„_)  - (m  - eJ/2,  and 

^ * 


or 


M-^x 


2 


1 

2 


As  (€j^  - m)/kT  grows  large,  (c  ) - kT,  and 

^ € 

X 


or 


rate  of  change  of  emitted  current  density  with  total  kinetic 
A/cm^/eV 


= Pn(£^)v^e 


energy, 
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These  functions  of  total  kinetic  energy  PN(ET)  and  J(ET)/DET  are 
distributions  of  electrons  emitted  because  of  positive  x-directed  velocity 
components  only,  not  because  of  total  velocities. 


Last  Two  Rows 


Under  the  columns  lie  two  rows  of  averaged  results  obtained  as  summed  approxima- 
tions of  equation  (4).  Although  the  emitter  electron  number  density  (NEM)  appears  first, 
the  remaining  properties  characterize  the  emitted  (or  about  to  be  emitted)  electrons. 


NEM 

NEE 

VXAV 

KEXAV 

KEXFL 

J 

KETAV 

KETFL 


TZERO 

TD 


electron  number  density  within  emitter,  complete  integration  of  eq.  (3), 

Q 

nem>  electrons/ cm 

number  density  of  the  emitted  electrons,  the  effective  integration  of 

O 

product  eqs.  (2)  and  (3)  over  all  v^^,  n^^,  electrons/ cm 

average  velocity  of  emitted  electrons,  eq.  (4)  with  u(v^)  = v^^,  (v^) , 
cm/sec 

average  kinetic  energy  based  on  x-directed  velocity  components  of 
emitted  electrons,  eq.  (4)  with  = mv^2,  {e^)  = m(  v^)/2,  eV/ 

electron 


average  flow  of  kinetic  energy  based  on  x-directed  velocity  components 
of  emitted  electrons,  eq.  (4)  with  u(Vj^)  = mv^2,  = m(v^)/2, 

(eV)(cm)/  sec/electron 


current  density  of  emitted  electrons,  j = ng^e^v^).  A/ cm 


average  total  kinetic  energy  of  emitted  electrons,  eq.  (4)  with 
u(v  ) = + <€  ) , + <€  ) ),  eV/electron 


average  flow  of  total  kinetic  energy  of  emitted  electrons,  eq.  (4)  with 
= Vs  < Sz)  < ^xS ) = < + '"x  < ^ ’ (eV)(cm)/ 


sec/ electron 

total  temperature  of  emitted  electron  stream,  T^  = 2(€rp)/3K,  Si 


static  temperature  of  emitted  electrons  radomized  about  drift  velocity 
<S>>  Tj5  = 2((€,j.)-m<v^)V2)/3/c, 
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Plots 


The  machine  plotted  the  last  four  columns  of  each  tabulated  set,  PN(EX)  against 
EPSX,  J(EX)/DEX  against  EPSX,  PN(ET)  against  EPST,  and  J(ET)/DET  against  EPST, 
and  tabulated  the  computed  results  on  each  figure. 

Figure  Index 
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FORTRAN  STATEMENT  FOR  CALCULATIONS  OF  THERMAL,  FIELD  EMISSION  WITH  A TERMINATED  IMAGE  POTENTIAL 

by  Susan  L Button  and  Annie  J.  Easley 
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o r>  o r»  r>  o n o r> 


MAINJM 

EXTERNAL  FORMULA  NfJMRCR  - SOURCE  STATFMFNT  - INTERNAL  FOPMtJLA  NUMBER  (S> 
r>ELTVX(N4*l)  = VX(N4.1)  - VX(N) 

MEAN  VALUE  OF  EPSLON  AND  VX 

fiy 

FPSLNM(N  + n = FPSLONfN^-n  - •St'AINCRINI  ,7B 

VXM(N)  * SQRT(2.0*EE>SLNM(N)/AC0N) 

COMPUTE  M/7  * VXM  ♦♦B 

,79 

VXMCU(N)=AC0N  /(2,0)  ♦ VXM(N)*»B  , BO 

70  CONTINUE 

RVX  = NVX 

SUMVX  ^ INDIVIDUAL  SUMS  , SUMVXI=  TOTAL  SUM  TO  INFINITY 

,Bl  ,B7 

SUMVXI  = 0.0  ,BB 

no  150  N * 1,300  ,84 

IF  (T(m  120,120,124  , R5 

120  IF(EPSL0N(N1-AMU(  jn  172,122, 167  , 86 

122  RVX(N)*I (4.0*3.1415  9*AC0N**?*B.61 727F-5+DEL T VX ( N ) ♦ t (AMJ( J)-FPSL  NM( 

IN)  )/8.61  727F-5M /(PARTI ♦♦?)) /PARTI  , 87 

GO  TO  140  ,88 

124  pART=4.0^3.14159*AC0N**2»8.61 727E-S*DELTVX (N) *T ( I ) , B« 

CONSTA  = (AMUU)  - E>SLNM(N))  / (8.61727F-5  ♦ T(D)  ,90 

IF  (CONSTA)  126,126,130  ,91 

126  IF  (ABS(C0NSTA)-88.028)  132,132,178  ,92 

128  RVX(N)  * 0.0  ,9R 

GO  TO  140  ,94 

130  IF(C0NSTA-88.028)  137,132,134  ,95 

132  PART?  = OLOGI 1.0DO+OFXP(CONSTA) ) ,96 

GO  TO  136  ,97 

134  PART?  * CONSTA  ,98 

136  IF  (PART?)  139,138,139  ,99 

138  PART?  = OEXP(CONSTA)  ,100 

139  RVX(N)  = PART*PART7/(PARTl**2) /PARTI  ,101 

140  IF  (N-n  142,142,144  , 107 

142  SUMVXIl)  PVX(l)  ,103 

GO  TO  146  ,104 

144  SUMVX(N)  * SUMVX(N-l)  ^ RVX(N)  ,105 

146  SUMVXI  = SUMVXI  ♦ RVX(N)  ,106 

NN  » N ,107 

150  CONTINUE 

CALCULATE  EPST 

,108  ,109 

162  IF  (Tim  200,350,700  , 110 

200  00  330  N = 1,NN  ,111 

SUMT  = 0.0  ,112 

SUMN  » 0.0  ,113 

JK  = I ,114 

IF  (FPSLNM(N)-AMU(J))  210,210,270  ,115 

210  XMIN  » 0.0  ,116 

XMAX  s 88.000  ♦ AKAPPA  ♦ T(l)  ^ AMUU)  - EPSLNM(N)  ,117 

CONST  UN  )=  -(EPSLNM(N)-AMU(  J)  ) / ( AKAPPA^Td  ) ) ,1)8 

IF  (CONSTKN )- 88.028)  7 30,230,220  , 119 

220  CIREYZ(N)  = AMU( J) -EPSLNM ( N) -88.  000* AK APPA#T (I ) , 120 
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non 


I. 


* 


* 


mainjm 

EXTERNAL  FORMULA  NUMBER 


SOURCE  STATEMENT 


INTERNAL  FORMULA  NUMBFR(S) 


XMIN2=CIRFYZ(N) 

XMAX2*XMAX 


KJJ  = 0 

rENTE(N)  = SIMPSl(  XMI N2 , XMAX 2 , A T NT, KJ J ) 

EYZ  (N)=CIRFYZ(N)-(CIREYZ(N)*CIREVZ(N)  )/(2.0^(  AMJ(  J)-FPSLNH(  N)))4^ 
1 AKAPPA*T( n*CENTE(N)/(AMU( J)-EPSLNMIN)> 


GO  TO  260 
230  KJ  ==  0 

ENTF(N)  = SIMPS1(XMIN,XMAX,AINT,KJ) 

IF  (ABSICONSTKN)  ) - B8.02R)  2*^  0 ,2<fO , 

240  FYZ(N)  = ENTE(N)  f CONSTUN) 


250 


260 


270 


280 

200 


300 


310 

320 


33  0 


350 


360 

370 


GO  TO  260 

OEN(N)  = ALOGU.O  ^ EXP  ( CON  STl  t N)  )) 

EYZ(N)  = ENTE(N)  / OEN(N) 

EPST(N)  = EYZ(N)  EPSLNM(N) 

GO  TO  330 

AX(N)  = (FPSLNM(N)-AMUUI  ) / ( AKAPPA*T  U ) ) 
TFMtN)  = EXP(*flOAT( JK)^AX(N) ) /FLOAT ( JK*JK J 
OFNOM(N)=  EXP ( -FLOAT  UK)* AX (N) ) /FLOAT (JK) 

IF  (TEM(N))  280,320,280 

IF  (Moot JK,2) ) 200,300,200 

SUMT  = SUMT  ^ TFM(N) 

SUMN  = SUMN  + OENOMIN) 

GO  TO  310 

= SUMT  - TEM(N) 

= SUMN  - OENOMCN) 

JK  + I 
270 

(AKAPPA*Tnn  ♦(SUMT  / SIJMN) 
FYZ(N)  «■  FPSLNM(N) 


SUMT 
SUMN 
JK  = 

GO  TO 
EY7(N)  = 
EPST(N)  = 
CONTINUE 
GO  TO  370 
00  360  N = 
EPST(N)  = 

continue 

00  380  N = 


1 ,NN 
(AMU( J) 

1,NN 


+ EPSLNM(N))  / 2.0 


* 121 
,122 
»125 
, 124 

,125 
,126 
,127 
,128 
,129 
,130 
,131 
, 132 
,133 
, 134 
,135 
,136 
,137 
, 138 
,139 
, 140 
,141 
,142 
,143 
, 144 
, 145 
, 146 
,147 
, 148 
,149 
, 150  , 

,152 
, 153 
,1*54 
, 155 


TABULATE  INOIVIOUAL  SUMS  / TOTAL  SUM 

THREE (N )=SUMVX(N) /SUMVXI 

380  CONTINUF 
C 

r START  CALCULATION  OF  0(N)  AND  OM(N) 

r 

00  650  N = l,NN 

IF  (E(M)  ) 400,570,400 

400  A=EPSLNM(N) 

410  BFTA  =( AMU( J )*PHI (K)-A) 

IF  (EVMAX  - EPSLNM(N))  560,420,420 
420  C0NST  = SQRT( BETA*? .0*Ar  0N) *2. 0*? , 0*3. 1 4 1 59* I . 6 0 2F-1 9/6 . 6242E -34 
CON  = ( BFTA/(  2 .0*F  ( M)  ) - ( FSORO  / (B*n*t>HT  { K)  ) ) ) 

ARO0T  = ( 4.0*F  SMALL/ f FSORO*F (M) -4,0 ♦PHI ( K ) ♦BET  A ) ) »*2 

8KuOT=PHnK  )*F  (M)*(  AMUU1-A) 

qnnT=ARrinT*RRnnT 

IF  (EVMAX-FPSLNM(N))  560,430,440 

430  XUN)=rnN 


tl57 
, 158 


,159 
,161 
, 162 
, 163 
, 164 
, 165 
, 166 
, 167 
,16B 
, 169 
,170 
,171 
, 172 


151 

rl56 

,160 


19 


r>  o n o o n 


MATNJM 

FXTERNAL  FORMULA  NUMBER  - SOURCE  SfAfEMENT  - TNfFRNAL  FORMULA  NUMBER ( S ) 

X?(N)*CON  , 

GO  TO  450  ,174 

440  Xl(N)=CON*n.O-SQRTn.O<-ROnT))  ,175 

X2(N)  =CON^(  l.04-SQRTn.O«'RnnT)  ) ,176 

450  TF(Xl(Nn  460,470,470  , t77 

460  XUNI=0.0  ,178 

470  ANTVL*(X?(N)-XUNn /lOO.  ,179 

X(l)  = XKN)  ,180 

DO  510  L = 1,101  , iRi 

CHFCK*U0-C6(  M)*X(L)  /BETA) -ESORO/(  8ETA»(4.0*X  (L)  i-ES  ORO/ PH  I ( < 1 M ,182 

IFCCHECK)  480,400,490  ,180 

480  ARGIL ) * 0*0  ,184 

GO  TO  500  ,185 

490  ARG(L 1*S0RT(CHECK)  ,186 

500  X(Lfl)  * X(L)  ♦ANTVL  ,187 

IF  rX(U  -X2(Nn  510,510,520  ,188 

510  CONTINUE  ,189  ,190 

520  CAIL  FNTGRL( 101, ANTVL,ARG, ANSI  ,191 

ANS(10n*ANS(l0n  , 197 

CONAN sABS(CONST*ANSI 1 01)  ) 

CALCULATE  OM(N)  USING  MEAN  VALUES 

♦ 193 

530  !F(CnNAN-88*028)  550,550,540  ,194 

540  OM(N)*0.0  ,195 

GO  TO  600  ,19^ 

550  OM(N)»EXP(-CnNST*ANS(10in  ,197 

GO  TO  600  ,198 

560  0M(N)«1,0  ,199 

GO  TO  600  , 700 

570  IF  (FPSL0N(N)-(AMU(JM“PHMK)  n 580,580,593  ,201 

580  DM(N)  » 0.0  , 202 

GO  TO  600  ,203 

590  0M(N)  » 1,0  ,204 

600  CONTINUE  , 205 

650  CONTINUE 


SUMOV  = INOVIOUAL  SUMS  OF  0( VXM)NCVXM)OELTAVX 
SUMOVl*  TOTAL  SUM  OF  0 ( VXM) N ( VXM) OE LTAVX I * NEF 
SUMTFN  » VXAV 
SUMFL  * KEXAV 
SUMTW  a KEXAL 


,206 

,207 

SUMOVm  = 0.0 

, 208 

SUMOV!  * 0.0 

,209 

00  750  N a l,NN 

,210 

ONVXM(N)=OM(N)*RVX(N) 

♦ 211 

IP  !N-n  700,700,710 

♦ 212 

700 

SUMOVU)  a ONVXMd) 

, 213 

GO  TO  720 

, 214 

710 

SUMOVCN)  a SUMnV(N-l)  ♦ ONVXMrN) 

f 215 

720 

SUMOVl  a SUMOVI  f ONVXM(N) 

, 216 

750 

CONTINUE 

,217 

,21  8 

00  800  N a 1,NN 

,219 

SfX(N)=SUMOV(N) /SUMOVI 

, 220 

800 

CONTINUE 

• 221 

,222 

n r»  n 


MAINJM 

EXTFRN4L  FORMULA  NUHRFR  - SOURCF  STAfFMFNT 


INTFR^AL  FORMULA  NUMBER ( S ) 


SUMTFN=0.0 
no  8BO  N * l,NN 

TENfNI  = ONVXMtN)  ♦ VXM(N)  / SUMOVT 
SIJM1FN=SUMTFN+TEN(N| 

8*50  CONTINUE 
SUMEL*0.0 
00  900  N * l.NN 

FL(N)  » ONVXM(N)  ♦ ACON  ♦ VXM < N I ♦♦2/2. 0/ S JMOVI 
SUMEL=SUMEL^EL<N) 

900  CONTINUE 
$UMTW=0.0 
on  9S0  N = I#NN 

TWFL(N)  - ONVXMIN)  ♦ ACON  ♦ VXMI N) ♦♦3/2. O/SUMOV T 
SIJMTW  = SUMTW^TWFLI  N) 

950  CONTINUE 


»?23 

,27^ 

»225 

,226 

.227  f 

,229 

,230 

,231 

♦ 232 

f233  , 

»235 

,236 

,237 

,238 


KETAV  , KETFL 

KETAV  = 0 

00  1000  N = 2,NN 

SKFTAVCN)  = EPSTINl  ♦ ( ST  X ( N ) -S I X < N-1 M 
KETAV  = KETAV  ♦ SKFTAVfN) 

1000  CONTINUE 

KFTFL  = 0 

on  1050  N s 2,NN 

SKETFUN)=  VXMIN)^EPST(N)  ♦(  SIX(N)-SIXfN-in 
KFTFL  = KETFL  ♦ SKETFL(N) 

1050  CONTINUE 
C 

C THIRT  = J 

C TZERO  , TO  , NEM 

C 

THIRT  = SUM0VI*SUMTEN^1.602E-19 
T7FR3  ?,0  ♦ KETAV  / (3.0  ♦ AKAPPA) 

TO  = 2.0  ♦ (KFTAV  -ACON^  SUMTEN**2/2. 01  / (3.0  ♦ AKAPPAl 

OONE=SUMOVI*(  2.0*8.61727F-5^T(n^SUMTFN/3.0  ♦ SJMTW) 
OTWn=OONE^1.602E-l9/4. 186 
NEM  :=  2.0  ♦ SUMVXI 

r 

C PNEX  , JEXn 

C PNFT  , JETO 

r 

on  1200  N = 1 ,NN 

CONST  1 (N  )*-( EPSLNM(N) -AMU( J)1 /(AKAPPA«T(  n ) 

PNEXC  N1*0NVXM(N1 /( ACON*VXM ( N) ♦OELTV X (N ) 1 
JEXDf  N)  = PNEX(N)6VXM (N)*1 . 6 02F-19 
IF  (T(I))  1130,1130,1100 

1100  IF  (CnNSTKN))  1110,1110,1120 

1110  IF(A8S(CON$Tl(Nn-17.000)  1 1 60 , 1 1 , 1 1 50 

1 1 20  IF(C3NST1 (N)-88.028)  1 160,1160,1130 

1130  PNET(N)  = 2.0  ♦PNFX(N1 
JETO(N)  = 2.0  ♦JFXO(N) 
on  TO  1700 

1160  0EN(N1  = OLOG ( 1. Ono+OEXP (CONSTl (N1  n 
S AVF( N)=( I .O^EXP( -CONSTl (Nil) 

PNFT(N)  = (PNEX( N )^AKAPPA*Tf 1 l^SAVF (Nl^nFN(N)  1 /FYZ( N1 


,239  , 

,261 
,262 
»263 
,266 
,265  , 

f 267 
,268 
,269 
,250 


t25l 
,253 
,256 
• ?*55 
,256 
t257 


,258 
f 259 
,260 
,261 
,262 
,263 
, 266 
, 265 
,266 
, 267 
,268 
,269 
,270 
,271 
, 272 


228 


236 


260 


266 


2 52 
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HAINJM 

FXTFRNAL  FORMULA  NUM8FR  - SOORCF  STATFMENT  - INTERMAL  Formula  N'JMBFR(«:» 
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rv.’  rvj  fv  cv 


<v  fvj  CM  r I rsi  {^  rvi  fv  f\j  *\  r.  fv  rv  <vi  (M  (V  fvj  c\  r t*'  rr  rr.  r cr  r'  ro  ^ rr,  ^ P"  rr  JJi;  pr,  rr,  tr  rr\ 


c C O C C C C C C c 

I— ‘ip-i 

c~:  r-'.  c>  c c~  r r rv  ^ p- 

rv'  rv  f\)  M rv  r.'  ro  rv  »\j 

CC  COCC'CO  cc 

Orr  (<‘l^f^J^c^aoc  rr^ 

o rA.  c I-'  rv.  c c o ^ ^ 

^(v  f\i  fv.  fv  rv  fv  f\,  cviPM 

► % »» 

c c c c c c c c c c 

Cp-  rsjtrir.  caoD  crj 

cc  M--CCCO 

<Ni  rv  fVJ  fv  fv  rvi  fv.  rv.  fv  <v 


C C C C 

»f  fVj  cc  c 
^ w r 
ro  (N.  r.'  (\ 

c c c c 

rsj  fv.'  p-  a 
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SIBFTC  PLHTl  DEBUG 

SUBROUTINE  PLOTl 

COMMON  /HP/  PNEXf  B«50)  , JEXOIOSD)  t EPS  LON  (350)  , EPS  T ( 350  ) , NN,  NCOUNT , 
1 NST0P,PNET(3501 ,JFTD(350) 

COMMON  N,AMUf 10) ♦ J,T( 1 0) ,I ,AKAPP& 

COMMON  FPSLNM(350) 

DIMENSION  P(IO)  t YSTOPOOO) 

DIMENSION  Y(500),X(500)  *YA(500)  ^XM500)  , YB  (500  ) ♦ XB{  500  ) , YC(  500  ) , 
1 XC(500) 


RFAU  JFXO  ,JETD 

WRITE  (6, 10) 

, 1 

. 2 

10 

FORMAT  (2HPT) 

DO  5 K*  NCOUNT,NSTOP 

t3 

Y(K)  = PNFXfK) 

,4 

X(K  ) = EPSLNM(K) 

♦ 5 

5 

CONTINUE 

f 6 

. 7 

KODE*BA 

, B 

NP  * NSTOP  - NCOUNT  ♦! 

f P 

P(  l)*NP 

.10 

CALL  SCALF  ( NP,Y( NCOUNT) ,KRSTR) 

♦ 1 1 

CALL  PLOTXY  (X(NCOUNT),Y  ( NC HUNT ) ,KOOE , P ) 

♦ 12 

WRITE  (6tll) 

. 13 

. ) 4 

1 1 

FORMAT  (2HPL,56X,17HPN(EX)  VS.  FPSLDN) 

WRTTE(6,10) 

,15 

.16 

on  15  K=  NCnUNT, NSTOP 

.17 

YA(K)*PNFT(< ) 

.15 

XA(K)«EPST(K) 

.10 

15 

CONTINUE 

.2P 

.21 

CALL  SCALE  (NP» YA(NCOUNT) »KRSTP ) 

.27 

CALL  PLOTXY  ( XA ( NCOUNT) , YA  ( NC OUNT ) ,KODE , P ) 

.20 

WRITF(6*21) 

.24 

.25 

21 

FOPMAT  (?HPL*57X»15HPN(ET)  VS.  EPST) 

WRTTF(6,10) 

.26 

.27 

00  25  K*  NCOUNT, NSTOP 

.?B 

YB(K)=JEXO(<) 

.20 

XB(K)=EPSLNM(K) 

, 30 

25 

CONTINUE 

. 31 

,3'? 

CALL  SCALF  ( NP ,Y8 ( NCOUNT ), KRSTR ) 

.30 

CALL  PLOTXY! XB(NCnUNT) ,YR  ( NC OUNT ) , KOBE , p ) 

. 34 

WRITE  (6,31) 

.35 

.36 

31 

FORMAT  (2HPL,53X,?2HJ(FX)/PELEX  VS.  FPSLON) 

WRITF(6, 10) 

.37 

.3« 

DO  35  NCOUNT, NSTOP 

. 30 

YC(K)*JFTD(K ) 

.40 

XC(K)»EPST(X) 

. 41 

35 

CONTINUE 

,4? 

. 4 3 

CALL  SCALF  ( NP , YC ( NCOUNT) , KR STP ) 

.44 

CALL  PLOTXY  ( XC ( NCOUNT ), YC  ( NCOUNT ), KODE , P ) 

.45 

WRITE  (6,41) 

.^6 

,47 

Al 

FORMAT  ( 2HPL , 54X,20HJ(FT) /DELFT  VS.  FPST) 

RFTIJRN 

• 4B 

END 

,49 
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MEMORY  MAP 


SYSTEM 

FILE  BLOCK  ORIGIN 

NUMBER  OF  FILES  • 2 

1.  UNIT05 

2.  UNIT06 

FILE  LIST  ORIGIN 
PRE-EXECUTION  INITIALIZATION 
CALL  ON  OBJECT  PROGRAM 
OBJECT  PROGRAM 


1. 

DECK 

•MAINJM* 

4 

2* 

DECK 

•AINT  • 

4 

3. 

DECK 

•PLOTl  • 

4 

4« 

SUBR 

•.IBSYS* 

5. 

SUBR 

••LRCON* 

6. 

SUBR 

•.lOEX  • 

7. 

SUBR 

•.JBCON* 

6. 

SUBR 

•.LXCON* 

9. 

SUBR 

•.lOOEF* 

10. 

SUBR 

••lOCSF* 

11. 

SUBR 

•.LXSL  • 

12. 

SUBR 

••FPTRP* 

13. 

SUBR 

•.ERAS.* 

14* 

SUBR 

•FOUT  • 

15. 

SUBR 

•FCNV  • 

16. 

SUBR 

•FIOS  • 

17. 

SUBR 

•FIOH  • 

18. 

SUBR 

•FSEL  • 

19. 

SUBR 

•FWRO  • 

20, 

SUBR 

•FROO  • 

21. 

SUBR 

•UN05  • 

22. 

SUBR 

•UN06  • 

23. 

SUBR 

•FLOG  • 

24. 

SUBR 

•FXPF  • 

25. 

SUBR 

•FSQR  • 

26. 

SUBR 

•FXP2  • 

27* 

SUBR 

•FXP3  • 

2B. 

SUBR 

•.IDE.  • 

29. 

SUBR 

•.I0E56* 

30. 

SUBR 

*.LOGER* 

31. 

SUBR 

•.NGOEF* 

4 

32. 

SUBR 

•.RWOOE* 

33. 

SUBR 

•.SQRTN* 

4 

34. 

SUBR 

•.TLEXP* 

35. 

SUBR 

•.XEXP  • 

4 

36. 

SUBR 

•.ZUOEF* 

4 

00000  THRU 
02720 


02750 

027  54 

02777 

03004  THRU 

03004 

41176 

41260 

00000 

00000 

00702 

02652 

52171 

52525 

52744 

55237 

55376 

56042 

56046 

56114 

57636 

60021 

60764 

61051 

61100 

61132 

61133 
61136 
61253 
61355 
61453 
61546 
61635 
61656 
61672 
61677 
61703 
61723 
61730 
61733 
61737 


02717 


70753 
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37* 

SUBR 

•PISTUG* 

* 

61743 

38. 

SUBR 

•PLOTXY* 

♦ 

62530 

39. 

SUBR 

•KHAR  • 

* 

64333 

40. 

SUBR 

•SCALE  • 

♦ 

64375 

41. 

SUBR 

•FOXP  • 

4t 

64546 

42. 

SUBR 

■OLOGIO* 

* 

65013 

43. 

SUBR 

•FDXl  • 

♦ 

65261 

44. 

SUBR 

•SIMPS!* 

65450 

45. 

SUBR 

•FNTGRL* 

♦ 

70636 

46. 

SUBR 

•//  • 

77212 

(♦  - INSERTIONS  OR  DELETIONS  MADE  IN  THIS  DECK) 

INPUT  - OUTPUT  BUFFERS  70754  THRU 


UNUSED  CORE 


77210  THRU 


77207 

77211 


BEGIN  EXECUTION 


REFERENCES 


1.  Morris,  J.  F. : Thermal,  Field  Emission  with  a Terminated  Image  Potential,  NASA 

TND-2784,  1965. 

2.  Dyke,  W.  P. ; and  Dolan,  W.W. : Field  Emission.  Vol.  vni  of  Advances  in  Elec- 

tronics and  Electron  Physics,  L.  Marton,  ed. , Academic  Press,  Inc. , 1956, 
pp.  89-185. 

3.  Bohm,  D. : Quantum  Theory.  Prentice-Hall,  Inc. , 1951. 

4.  Sommerfeld,  A. ; und  Bethe,  H. ; Elektronentheorie  der  Matelle.  Handbuch  der 

Physik,  vol.  24,  pt.  2,  Springer-Verlag  (Berlin),  1933,  pp.  333-622. 

5.  Avak’yants,  G.M. : Theory  of  Electron  Emission  from  a Metal  in  an  Electric  Field. 

Trudy  Fiz. -Tekh.  Inst.  Akad,  Nauk  Uzbek  S.S.R. , vol.  6,  1955,  pp.  43-53. 
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T = 0.  E » 0.1C000002E  08  PHI  * 4.00  AMU  * 10.00  EVHAX  * 12.8901 

EPSX  EPST(VXH)  XI  X2  P(VXM)  PNOV  PNIEXJ  JlEXl/OEX  PNIETI  J(FT|/0ET 

9.4C000E  00  9.67500E  00  0.  4.57280E-07  2.17006E-28  2.36213E-05  7.86090E-08  2.28377E-10  1.57218E-07  4.5S754E-18 

9.5000CE  00  9.72500E  00  0.  4.47108E-07  2.03995E-27  l.88066E-04  6.21957E-07  1.81656E-17  1.24391E-06  3.63313E-17 

9.52000E  00  9.7550CE  00  0.  4.41001E-07  7.74361E-27  2.98940E-04  2.09fc74E-06  6.14339E-17  4.19347E-06  l.22868E-lf. 

9.5400CE  00  S.765C0E  00  0.  4.38965E-07  1.20581E-26  4.64371E-04  3.12843E-06  9.17586E-17  6.25686E-06  1.83917E-16 

9.560CCE  00  9.77500E  00  0.  4.36928E-07  1.87600E-26  7. 1053AE-04  4.85520E-06  1.366B3E-16  9.31041E-06  2.73366F-18 

9.5E000E  00  9.78500E  CO  0.  4.34891E-07  2.91664E-26  1.07586E-03  6.908A2E-06  2.03058E-16  1.381 72E-05  4.06117F-16 

j 9.600C0E  00  9.79500E  00  0.  4.32854E-07  4.53007E-26  1.61632E-03  1. 022068-05  3.00717E-16  2.04411E-05  6.01434E-16 

t 9.6200CE  00  9.80S00E  00  0.  4.30816E-07  7.02943F-26  2.41322E-03  1.50702E-05  4.43B68E-16  3.01404E-05  8.87737E-16 

9.640CCE  00  9.81500E  00  0.  4.28778E-07  1.09009E-25  3.58442E-03  2.21485E-05  6.53029E-16  4.42971E-05  1.30606E-15 

9.66000E  00  9.82300E  CO  0.  4.26739E-07  1.68856E-25  5. 298796-03  3.24202E-05  9,568726-16  6.484C5F-05  1.91374E-15 

I 9.680006  00  9.835C0E  00  0.  4.24700E-07  2.61344E-25  7.797986-03  4.726186-05  1.39636E-15  9.45236F-05  2.79272F-15 

9.7000CE  00  9.8450GE  00  0.  4. 226616-07  4.041806-25  1.142516-02  6. 859176-05  2. 028656-15  1.371836-04  4. 057316-15 

I 9.72000E  00  9.8550CE  00  0.  4.206226-07  6.24529E-25  1.66625E-02  9.90462E-05  2. 932396-15  1.980926-04  5.364786-15 

I 9.740006  00  9.065006  00  0.  4.185826-07  9.642156-25  2.418346-02  1.422266-04  4.215116-15  2.044526-04  8.430236-15 

9.7600CE  00  9.875006  00  0.  4. 165416-07  1.487486-24  3. 491536-02  2.029486-04  6.020916-15  4.058976-04  1.204186-14 

, 9.7800C6  00  9.885006  00  0.  4.145016-07  2.292746-24  5.011796-02  2.074976-04  8.537966-15  5.749946-09  1.707596-14 

! 9.800C06  00  9.895006  00  0.  4.124596-07  3.530756-24  7.147206-02  4.038266-04  1.20049E-14  8.076526-04  2.400996-14 

9.820006  00  9.905006  00  0.  4.104186-07  5.433756-24  1.011756-01  5.617186-04  1.671586-14  1.123446-03  3.343166-14 

9.8400CE  00  9.915006  00  0.  4.083766-07  0.355446-24  1.420016-01  7.720456-04  2.299826-14  1.544096-03  4.599636-14 

9.860CC6  00  9.925006  00  0.  4.063336-07  1.283826-23  1.972936-01  1,045646-03  3.117986-14  2.091276-03  6.235966-14 

I 9.860CC6  00  9.935006  GO  0.  4.042916-07  1.971316-23  2.708016-01  1.390106-03  4.149326-14  2.780196-03  8.298656-14 

9.900006  00  9,945006  00  0.  4.022476-07  3.025146-23  3.661546-01  1.803216-03  5.307096-14  3.606426-03  1.077586-13 

9.920006  00  9.955CCE  00  0.  4.002046-07  4.638866-23  4.856666-01  2.260096-03  6.759856-14  4.520186-03  1.351976-13 

9.940C0E  00  9.965006  00  0.  3.981596-07  7.111826-23  6.280346-01  2.692276-03  0.060606-14  5.38454F-03  1. 612126-13 

9.960006  00  9.975006  CO  C.  3.961156-07  1.089996-22  7.837336-01  2.944446-03  8. 024466-14  5.88880F-O3  1.76489F-13 

9.980006  00  9.985006  CO  0.  3.940706-07  1.670826-22  9.267996-01  2.705516-03  0.116556-14  5. 411036-03  1.62331F-13 

l.COOOOe  01  9.995006  CO  0.  3.920246-07  2.566606-22  1.000006-00  1.384296-03  4.157056-14  2.760596-03  0.314096-14 

N6M  = 0.143657326  24  NEE  » 0.378218166-03  VXAV  * 0.186675376  09  K6XAV  » 0.99C730836  01  K6XFL  * 0.184949146  10 

J » 0.1131086-13  K6TAV*  0.9953656  01  KETFL*  0.1858126  10  TZ6R0  * 0.7700556  05  TO  » 0.3594036  03 

I 


T X 0.  6 a 0.100000026  08  PHI  » 4.00  AMU  = 10.00  EVMAX  * 12.8901 

NEM  a 0.143657326  24  N66  * 0.378218166-03  VXAV  = 0.186675376  09  K6XAV  * 0.990730036  01  KEXFL  * 0.184949146  10 

J » 0.1131086-13  KETAV*  0. 995365E  01  KETFL=  0.185812E  10  TZERO  * 0.770055E  05  TO  * 0.3594036  03 

Figure  2.  - Sample  of  complete  IBM  output  sheets. 
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T * 
NEH  « 
J = 


0.3C0000C0E  0^ 
0.49014<»90E  22 
O.A7l3^6E  10 


E » 
NEE  > 
KETAV* 


0.1C000002E  09 
0.84568L50E  21 
0.897721E  00 

Figure  3.  • 


PHI  = 2.00  AHU  * UOO  EVMAX  = 1,0053 

VXAV  * 0.3A791286E  08  KEXAV  0 . <»  33  7033  IE  - 00 

KETFL*  0.361517E  08  TZERO  = 0.69A513E  0<<  TO  = 

• Characteristics  of  thermal,  field  emission. 


KEXFL  * 0.21864102E  08 

0.<»28283E  0^ 


I 

I 


I 


I 
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T ■ 0.30000000E  0^  6 ■ C.ICOOOOOZE  09  PHI  • 2.00  AMU  • 5.00  EVMAX  » 5.0053 

NEM  - 0.5092984JE  23  NEE  • 0.369082UE  21  VXAV  - 0.U9A4921E  09  KEXAV  » 0.410A9357E  01  KEXFL  • 0.50097755E  09 

J - 0.706267E  10  KETAV«  0.A71579E  01  KETFL-  0.569826E  09  T2ER0  ■ 0.36A833E  05  TO  • 0.510114E  OA 

Figure  3.  - Continued. 


32 


T - 

0.300000COE  04 

6 » 

C.1C000002E 

09  PHI 

NEM  • 

0.14377540E  24 

NEE  « 

0.24587175E 

21  VXAV 

J » 

0.706308E  10 

KETAV-  0. 

974862E  01 

KETFL*  0 

2.00  AHU  s 10,00 

0.17931756E  09  KEXAV 
198E  10  TZERO 

Figure  3.  - Continued. 


eVMAX  » 10.0053 

• 0.91592062E  01 


» 0.754192E  05 


KEXFL  « 0.16^S5962E  10  ' 

0.469611E  04 
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T . o.soooooooe  04  e . c.iooooooze  o«  phi  . z.oo  *hu  - 15.00  evh»x  . 15.0051 

NiH  ■ 0.26404021E  24  NEE  » 0.19754E24E  21  VXAV  ■ 0.221le09TE  09  XEXAV  • 0.t4171910E  02  XEXPl  • 0.1I677492E  10 

J - 0.7063C4E  10  KEt«V>  0.147565E  02  KE7FL«  0.329644E  10  72ER0  • 0.1141&2E  06  TO  ■ 0.460799E  04 
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Figure  3.  - Continued.  ^ 
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Figure  3.  - Concluded. 
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aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl- 
edge of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof.*^ 

— National  Aeronautics  and  Space  Act  of  1958 
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